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Abstract 

In the interaction of an acoustic field with a moving airframe we encounter a canonical 
initial value problem for an acoustic field induced by an unsteady source distribution, <?(<,*) 
with q = 0 for t < 0, in a medium moving with a uniform unsteady velocity U(t)i in the 
coordinate system x fixed on the airframe. Signals issued from a source point S in the 
domain of dependence V of an observation point P at time t will arrive at point P more 
than once corresponding to different retarded times, r in the interval [0, t]. The number of 
arrivals is called the multiplicity of the point S. The multiplicity equals 1 if the velocity U 
remains subsonic and can be greater when U becomes supersonic. For an unsteady uniform 
flow U(t)i, rules are formulated for defining the smallest number of I subdomains V, of V 
with the union of Vi equal to V. Each subdomain has multiplicity 1 and a formula for the 
corresponding retarded time. The number of subdomains Vi with nonempty intersection is 
the multiplicity m of the intersection. The multiplicity is at most I. Examples demonstrating 
these rules are presented for media at accelerating and/or decelerating supersonic speed. 
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1. Introduction 


In this paper we study the propagation of acoustic waves induced by a spatial source 
distribution q(t, x) in a uniform medium moving with a time dependent velocity 17(f)* 
relative to the coordinate system x = (x, y, z). This study is an outgrowth of recent exper- 
imental and theoretical investigations of a model acoustic-panel interaction problem. See 
for example [1] and [2]. The acoustic source distribution could simulate the noise from jet 
exhausts, turbulent boundary layers and/or the airframe panel oscillation. For the analysis 
of the panel oscillation it is convenient to work in the coordinate system x = (x,y, z), sta- 
tionary with the airframe, while the flow field is moving with the uniform unsteady velocity 
In this system, the acoustic potential $(f , x, y, z) is governed by the convective wave 
equation with variable coefficients, the velocity 17(f) and the acceleration 17(f), 

+ Ud x +2Udf X + U 2 dl x — C 2 A]$(f,x,y,z) = q(t,x,y,z) , (1.1) 

where A denotes the Laplacian in (x,y, z) and C stands for the speed of sound of the 
medium at the ambient condition. For an initial value problem, we impose 

$ = 0 , d t $ = 0 , at f = 0 , (1.2) 


or 


9 = $ = $ ( = 0, for f < 0 . (1.3) 

In general, the solution $(f , x) is constructed indirectly via the solution of the correspond- 
ing acoustic field $ in the coordinate system x with the medium at rest [3]. These two 
coordinate systems, which coincide at f = 0, are related by translation, x = x + *A'(f), 


l. e.. 


x = x - X(t) , y = y , 2 = 2 


with X 


(() = f U(t') 

Jo 


dt' 


(1.4) 


From hereon we omit the bars over y and z for the coordinates of x. Figure 1 shows the 
coordinate x with the medium at rest moving with velocity 17(f)* relative to the coordinate 
x. We use / and / to denote the same quantity as functions of f, x,y,z and f,x,y, z 
respectively with t = f . Thus we have 

/(f,x,y,z) = f(t,x,y,z) with d-J = (d t + Ud x )f and V/ = V/, (1.5) 


where V and V denote the gradient operators in x and x. In particular, if / denotes the 
unsteady source distribution q in the coordinate system, x, the corresponding distribution 
q in x is given by (1.5). In the barred system, the acoustic potential $ is governed by the 
simple wave equation, 


[C 2 d%t- A] §(t,x,y,z) = q(t,x,y,z) , (1.6) 

where A denotes the Laplacian in (x,y, z). With homogeneous initial conditions (1.2) 
applied to the barred variables, the potential at point P(x) at instant f, is given by the 
explicit formula [3], 


$(t,x,y,z) = J j J_^d£dr]d( 


g(r, J , ,y , ,z') 
R 


(1.7) 


where ^,y,C denote the coordinates of a source point S(x') relative to point P(x), 
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(1.8) 


x' = x + { , y' = y + y , J = 2 + C , 
i? denotes the distance from 5 to P and r the retarded time, 

R = |x' — x| = \j i 2 + o 2 and r = t — RJC where a = \A ? 2 + C 2 • (1-9) 

We now identify the domain of integration V as the domain of dependence V of point P 
at time t. In the coordinate system, £, r), £, the point P is the origin of the system and the 
domain of dependence T> is a spherical ball of radius Ct centered at the origin, P, i. e., 

®={(f,>?,C)|0 <«<«}. (1.10) 

It is covered by a family of concentric spheres of radii R € [0, Ct]. The signal from a point 
S in V on a sphere of radius R arrives at the origin P along the radial line SP at time 
t when the signal was initiated at the retarded time t = t — R/C. Since the domain of 
dependence is axi-symmetric with respect to the £ axis, we show in Fig. 2 the cross section 
of the domain P, its boundary B and the family of concentric spheres in a meridian plane, 
the plane. We define the multiplicity of a point by the number of retarded times the 
signal from the point reaches point P at time t. Thus the multiplicity of a point in D 
or the multiplicity of T>, in the barred coordinate system with the medium at rest, is one 
and hence the domain of integration V for the integral on the right hand side of (1.7) is 
equal to the domain of dependence T>. Note that the above description of the domain of 
dependence (1.10) in the barred coordinate system is independent of the velocity of the 
medium relative to the coordinate system x. 

To identify the source distribution q at point S(x' ,y', z') at a retarded time r with 
the distribution q at the corresponding point S(x',y',z') in the moving frame, we relate 
the coordinate x' to x' by (1.4) at time r, i. e., 

x' = x' + X(r) and q(T,x' ,y\z') = q(r,x' ,y' ,z') . (1.11) 


By using (1.4) and (1.5), the left hand side of (1.7) is identified as $(t,x,y, z) and (1.7) 
becomes an implicit formula for $ at point P(x, y, z). Since the dependency on the velocity 
U(t)i of the medium appears only via the transformation (1.11) for the source distribution, 
the solution (1.7) of the initial value problem, (1.6) and (1.2), is valid for all Mach numbers 
including the transonic range. 

To express the potential $ as an integral in the coordinate system x, we change the 
integration variables ^,77, ( in (1.7) to £,77, £ which are the coordinates of point S relative 
to P. Using (1.4), (1.8) and (1.11), we have £ as function of £ and a or £ and r, 


£ = x' — x 


= £ + X(r) - X(t) = £ - 



where t = t — R/C , R = i/£ 2 + cr 2 and <7 = -y/77 2 + £ 2 . 

The Jacobian determinant of the coordinate transformation is 


d(£,?7, C) 

d(£,f?,C) 


5 £ 

% 


R — M£ 

R 


1 — M cosrp. 


( 1 . 12 ) 


(1.13) 
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Here R, 8 , ip with 0 < ip < it represent the spherical coordinates while £, a, 8 the cylindrical 
coordinates with 77 = a cos 8 and ( = cr sin0. Note that in the transformation (1.12), 77 and 
C and hence <7 and 8 remain the same. The cylindrical coordinates in the £, 77 , ( space are 
£, a and 8. 

The spheres S of constant r in a subsonic interval, where M(r) < 1 , cover a subsonic 
region in V. The spheres in a supersonic interval, where M(r) > 1 , cover a supersonic 
region. In a subsonic region, the Jacobian remains positive. In a supersonic region, there 
is a critical surface, J , on which the radial velocity component is sonic, 

J = 1 — M(r) cos = 0 or a = £\/ M 2 — 1 > 0 . (1-14) 

The transformation ( 1 . 12 ) maps the domain of dependence V of the origin P of the 
£, 77 , ( space to the domain of dependence T> of the origin P of the £, 77 , ( space. The critical 
surface J in a supersonic region is then mapped to surface J. Since the Jacobian J changes 
sign across 3 1 the neighborhoods on the two sides of 3 axe mapped to the same side of 3 • 
Thus the inverse transformation of (1.12) is not single valued and hence the multiplicity 
of a point in V can be greater than 1. For a given point (£,cr) and time f, (1.12) becomes 
an implicit equation for £ because the lower limit of the integral t — t — Rj C and R are 
related to £ by (1.9). Using (1.9) and (1.12), we get an implicit equation for R(t, 

i ? 2 = [£ + fuit'Wtf+o 2 . (1.15) 


The root(s) of (1.15) in [0, Ct) depend on the unsteady velocity U(t'). 

Using (1.11) to (1.13), we convert (1.7) to a volume integral in the £, 77 , ( space, 


§{t,x,y,z) = 


- Wl *** 

-him. 


g(u* + f,y + »?,z + C) 


d^drjd( 


\R-MZ\ 

g(^a + £,y + ?7,z + 0 
\Ri-MiZ t \ 


(1.16a) 

(1.166) 


To evaluate the integral over the domain of dependence V in (1.16a) we need to relate R 
and l to (, o and t , that is, to find all the roots, i?, of (1.15) in [0, Ct] or the retarded time 
Ti in [0,t] and identify the domain Vi in the 77 , ( space for each R{. Then $ is written as 
the sum of integrals over the smallest number I of domains V t in (1.16b). It is the purpose 
of this paper to formulate the rules identifying these domains Vi for a medium moving at 
an unsteady speed, U{t). First, we review briefly the solutions given in [2] and [3] for the 
following three special cases; 1 ) moving at an unsteady subsonic speed, 2 ) moving at a 
constant supersonic speed, and 3) moving with constant acceleration. 

1.1, Unsteady Subsonic Speed For this case, we have M(r) < 1 , r € [0, t], the forward 
and backward speeds of propagation C ±U are positive and hence the family of concentric 
spheres covering V is mapped to a family of nested spheres covering the spherical ball V 
bounded by Sq. Hence 1=1 and the domain of dependence V has multiplicity 1 . This 
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result reduces to the well known result for a steady subsonic flow, at a constant M < 1, 
for which (1.15) has a unique and explicit solution for £,? 7 ,£ in V and the denominator in 
(1.16) becomes the pseudo- distance [3], 

R - M£ = R* = v^ 2 + (1 — M 2 )< 7 2 . (1.17) 

The result that the domain of integration in (1.16) is equal to V remains valid even 
when there is a sonic interval, where M(r) = 1 for 0 < ti < t < t2 < t, during which, the 
nested spheres will be in contact at one point on the £-axis with £ = C(t— r) — X(t)+X(r) = 
C(t — ti) — X(t) + X(ti ). Although the multiplicity at the point of contact is infinite, the 
integral in (1 .16) remains finite because the acoustic field is initiated at t = 0. 

1.2, Constant Supersonic Speed For the case of a medium moving at constant super- 
sonic speed, M > 1, the solution is well known [3]. The domain of dependence V of point 
P lies in the backward Mach cone from P. The cone is 


£ = - \/(M 2 - 1)<t 2 . (1.18) 

To evaluate the integral we need to relate R and £ to £, a and t, The domain T> is covered 
by a one parameter family of spheres of radius R = C(t — r) centered on the £-axis at 
£ = —MR for r € [0, £]. The Mach cone is the envelope of these spheres. The circle of 
tangency of a sphere with the cone divides the sphere into a large and a small surface, 
convex and concave away from the cone respectively. Fig. 3 shows the cross section of 
the Mach cone and the spheres in the meridian plane, the £<r plane. As an initial value 
problem, we have r > 0 or R < Ct , thus the domain of dependence V is bounded by the 
backward Mach cone (1.18) and the convex part of the spherical surface So with radius 
Ct centered at £ = —Ut and azimuthal angle ip € [7t/2 — a, i r], where a is the Mach angle. 
Note that the mirror image of the domain of dependence V with respect to the plane £ = 0 
is the domain of influence of a source S at the origin for the duration [0, t] described in [3]. 
With a constant M > 1, (115) reduces to a quadratic equation for R and the two roots 
are [3], 

R ± = -[MZ±R*)/(M 2 -1) , (1.19) 

where R* = y / £ 2 — (M 2 — l)a 2 . The denominator of the integrand in (1.16) becomes R* 
and the integral can be written as the sum of two integrals over the domains V~ and V + . 
Equation (1.16) becomes, 


$((, x,y,z) = -S [J J J vt + J J J y _ 1 


qijxx + £,y + T?,z-K) 

R* 


d(drjd( 


( 1 . 20 ) 


with r = t — R± in V ± and 1 = 2. The domain of integration V + is bounded by the Mach 
cone and convex part of the sphere So and is equal to T>. The domain of integration V “ 
is bounded by the concave part of the sphere So with <p £ [0, a]. Note that the superscript 
of V gives the sign of J in V±. Since V~ € V + = T>, the multiplicity of T> fl V~ is two 
while the multiplicity of the complement, the spherical ball inside So, is one. The meaning 
of and the propagation of signals from point S at r* to P at time t are described in 
Appendix A. 
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Without making use of the explicit solution (1.19) of (1.15), the results in [3] quoted 
above were obtained by a different method in [2], It begins with the partition of the domain 
of dependence the spherical ball (1.10), into two nonoverlapping subdomains, V ± , by the 
critical surface J , (1.14), which is the circular cone, ip = arccos(l/M) = ir/2- a. Shown 
in Fig. 4 is the partition of V in a meridian plane by the critical curve J , the radial 
line PT, into two subdomains V~ and V + , with 0 < ip < 7r/2 — a and 7r/2 — a < ip < it 
respectively. The image of PJ in the £<7 plane in Fig. 3 is the backward Mach line PT and the 
images of V± are the domains of integration V± described in the preceding paragraph.. The 
procedure for the identification of the domains of integration in the £, 77, ( space, introduced 
in [2], demonstrates the important role played by the critical surface, which does not exist 
in an unsteady subsonic region. Since this procedure does not need the explicit solution 
of (1.15) or (1.12), it is suitable for generalization to an unsteady speed U(t) with several 
supersonic intervals. As an illustration or rather a clue to the generalization, we mention 
the special case of constant accelerating motion considered in [2]. 

1.3 Constant Accelerating Motion The interaction of an acoustic field with a baffled 
panel in the horizontal xz plane was studied in [2], for the panel moving with constant 
acceleration to a supersonic speed, i. e., U = Ct . With M = t for r 6 [0,t] and t > 1, 
there is a supersonic interval 1 < t < t. This special case shows that in addition to the 
partition of V by the critical surface j in a supersonic region, further partitions may be 
needed, related to the geometrical properties of J and its image J . See Figs. 4 and 
5. This case is discussed in §4 as one of the examples for the application of the rules of 
partition formulated in §3. 

In the next section, we study the geometry of the critical surface J and the corre- 
sponding singular points on the image PJ ■ We then use these studies to formulate in §3 the 
rules for the partition of the domain of dependence, V, into nonoverlapping subdomains 
Vj ' s for an unsteady velocity, U(r) , r € [0, t]. and show that the partition has the required 
property that the transformation (1.12) from a Vj to its image V } in the £,77, ( space is one 
to one. In particular, the value of R or r defined by (1.15) for a point in Vj is uniquely 
defined. Thus the domain of integration V in (1.16) is the sum of those images Vj' s. Since 
those images may overlap each other in the £, 77, ( space, signals from a point 5 in the space 
lying in ui of those images will reach point P 7n-times. Hence vn is the multiplicity of point 
S. The rules for counting the multiplicity of a point in D are also stated in §3. The anal- 
yses in §2 and 3 justify the assertions made in [2] for the case of a constantly accelerating 
speed. Additional examples with unsteady velocities, accelerating and/or decelerating at 
different ranges of Mach numbers are presented in §4 to demonstrate the application of 
the rules of partition of V and to show the singularities of J in the 77. ( space. 

2. The Geometry of the Critical Surface and its Image 

In § 1 , the domain of dependence V of point P at time t and the critical surface J in 
the coordinate system (£, 77, £) with the_ medium at rest are defined by (1.10) and (1.14). 
The transformation from the system (£ , T7,0_to (£,77,0 defined by (1.12) at time r is a 
pure translation of (X(t) - X(t)) along the £ axis. We reemphasize the following points) 

f Numbers in curly brackets refer to statement numbers in §2. 
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{1} The domain V is a spherical ball of radius Ct centered at point P and for each point 
on a concentric sphere of radius R € [0, Ct], the retarded time is t — R/C. Thus V 
has multiplicity one, independent of the source distribution and the velocity U(r)t of 
the medium. 

{2} The critical surface J , on which the Jacobian determinant J vanishes, appears in V 
only in a supersonic region. The surface is axisymmetric with respect to the £-axis, 
is independent of the source distribution but depends on the velocity U(t) of the 
medium. See for example Figs. 2 and 4 

{3} The transformation (1.12) relates the axial variable £ to £ while the other two variables 
77 and £ or a and 0 remain unchanged. The transformation is also independent of the 
source distribution but depends on U(t). Thus the images of the spherical ball T>, 
the concentric spheres of constant r and the critical surface J in the ( space are 
axisymmetric with respect to the £-axis. 

{4} A spherical cap of radius R centered at the origin of the £,?7,C space is mapped to 
the same spherical cap centered at i(X(r) — X(t )) in the £,77,^ space with the same 
retarded time r = t — R/C. See for example Figs. 3 and 5. 

Because of the axisymmetry, it suffices to study the geometry of the surfaces J and its 
image J in a meridian plane, the (<r plane, and the corresponding £<7 plane. We retain the 
same symbols V, J , Vj for their cross sections in a meridian plane and the same unbarred 
symbols for the images. The geometrical properties of the critical curve J and its image 
J will be described in § 2.1 and § 2.2 respectively and be employed in §3 to formulate 
the rules for the partition of V into nonoverlapping subdomains Vj, j = 1, 2, ■ • ■ such that 
the mapping of a Vj to its image Vj is one to one. 

2.1. The critical curve J in a meridian plane, the £<7 plane 

In the meridian plane, the two spherical coordinates R and ip become the polar coordi- 
nates, and the domain of dependence V in the plane is a semicircular disc with 0 < R < Ct 
and 0 < ip < 7r. Let us consider the signals initiated in a supersonic interval, i. e., 
M(r) >1, r € and M = 1 at if > 0 and at t* if t* < t. The signals come 

from points in the supersonic ring C(t — t*) < R < C(t — <*) in V. The critical curve J 
in the polar coordinates is defined by, 

cosip = 1/M(r) and R = C{t — r) , (2-1) 

and in the rectangular coordinates (V a) by 

£ = R cos ip = C(t — t)/M and <7 = Rsinip = C(t — r)[l — M -2 ] 1 / 2 , (2-2) 

with t 6 [<*,<*]. See for example Fig. 4, where = 1 and t* = t. Since the Mach angle 
a(r) = t/2 — ip we have from (2.1), 

{5} The radial line to point Q(R,ip) on J is orthogonal to the Mach line at Q. The 
polar angle ip of Q reaches a local maximum (minimum) when M or U reaches a local 
maximum (minimum). 
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The rates of change of a and £ are 


„ v RM' - CM(M 2 - 1) 7U . -CM - RM' 

J _ (r)= w-i ’ * (T) = — w — 

and the slope of J is 

do _ [RM 1 - CM(M 2 - 1)] 

~ -(CM + RM')y/AP - 1 ' 


(2.3a) 


(2.36) 


From hereon we use the prime to denote the r-derivative. From (2.3), we see that o'(t) < 0 
when M' < 0, therefore, o' and the slope vanish only when M 1 > 0 while £ ; (t) < 0. Thus 

{6} The ordinate a of J attains an extremum only in an accelerating subregion, while £ 
decreases. 

Since o remains unchanged in the transformation from the £ a plane to (o plane, the 
mapping of a neighborhood of the curve J to its image is defined by the mapping of a 
horizontal line segment at constant a. A small segment which crosses over J at a point 
with nonzero slope, or o' ^ 0, will cross over J only once. 

Now consider an extremum point, say K e , i. e., o' e = 0 at r e , we have 


(' e = -CM t and o" 


R e M" - 3 CM e 2 M' 

T 


(2.4a) 


We have assumed that U"(r) or M" exists, therefore, o" exists with M e > 1 . This is a 
realistic assumption and will simplify our discussion that follows but is not essential. 

When the extremum point is neither a maximum nor a minimum, o' does not change 
sign as r crosses r e and o " = 0. A horizontal segment close to point K e will intersect J 
only once. 

When o" ^ 0 at a critical point say K\, it is either a local maximum or minimum 
point. The curve J near K\ behaves as a parabola with (o — o\ )o" > 0. For a given <r, 
there are two roots, 


£ - £l = ±CM l \ 2{< ’ /’V 72 + 0(a - < 7 ,) . (2.46) 

a i 

Thus a horizontal segment near point K\ crosses over J twice if it is below (above) K\ 
where o is a local maximum (minimum). This point will be elaborated in {11a} and 
{116}. We need to locate the local maximum and minimum points on J. Let the n-th 
point and the corresponding time be denoted by K n and r n with r„ in decreasing order, 
i. e., r n > T n + 1 , for n — 1,2 •••,1V. The corresponding radii R n = C(t — r n ) are in 
increasing order. These N points divide J into N + 1 segments. The segment from point 
K n - 1 to K n , is denoted by J n with r 0 = r* and r^r+i = r*. 

If the supersonic interval ends at t* = t, we have R(t*) = 0. If t* < t, the interval 
ends with M(t*) = 1 due to the continuity of M(r). For either case, we see from (2.1), 
(2.2) and (2.3), that 
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{7a} For the supersonic region, the critical curve 3 ends at r = t* on the ^ axis with 
cr'(t*) <0. If J has one or more extremum points, N > 1, the first one K\ has to be 
a local maximum. Therefore, K n is a local maximum if n is odd and is a minimum if 
even. 

{76} If the supersonic interval begins at <» > 0, we have M{t *) = 1 and cr'(t) — ► +oo as 
t — + t* + 0. Hence K,\ is a local maximum and N is odd. If N is even, it is necessary 
that t* = 0. 

If the supersonic interval begins at f* = 0, we have M( 0) > 1 and the sign of cr'(0 + ) depends 
on the initial acceleration or M'( 0 + ). For example, in the case of constant supersonic speed, 
M'(t) = 0 and hence a' < 0 as shown in Figs. 2. On the other hand, the sign can be 
positive, if the interval begins with an acceleration at a low supersonic speed. See for 
example Fig. 4. 

Because of {7}, the number N of extrema of j will be even or odd depending on 
whether is negative or not. Thus 

{8a} If the ordinate of the critical curve J decreases at the beginning of a supersonic 
interval, i. e., a'(t, + 0) < 0, N has to be even. The curve has either no extremum or 
pairs of maxima and minima, with the last one being a local minimum apj at Kjf. 

{86} If the ordinate of the critical curve J is not decreasing at the beginning of a supersonic 
interval, i. e., cr'(t* + 0) > 0, N has to be odd. This is the case when f* > 0. The curve 
J has at least one local maximum a at point K\ and time Ti . If there are additional 
ones then they come in pairs of minima and maxima. The last one is a local maximum 
at Kp;. 

Note that {76} is included in {86}. Since J has an extremum a only when M' > 0, {8a} 
and {86} are applicable to an accelerating supersonic interval. 

From (2.4a) for a" at an extremum point and {6} that M' > 0, it is clear that cr" < 0 
when M" < 0 and the extremum point can not be a local minimum. Therefore 

{8c} In a supersonic region with constant or decreasing acceleration, there can be at most 
one local maximum and no local minimum. 

Note that for the case of a constant supersonic speed [3] reviewed briefly in §1.2, 
the critical curve j which is a radial line does not have a local extremum a. This result 
is in agreement with {8a} since a'(r) = -1/M ,r 6 [0,i]. For the case of constantly 
accelerating speed from subsonic to supersonic [2] 2 reviewed briefly in §1.3, we can use 
{86} and {8c} to conclude that the critical curve J has only one maximum point and no 
minimum as shown in Fig. 4. 

Figure 6 gives another demonstration of {8a}. It is the case of accelerating at a 
constant rate from a low supersonic speed, with M(r) = 1.15 + r 2 for 0 < r < t = 2.25. 
The critical curve j begins at point f on the semicircle B with M( 0) = 1.15 and <r(0) = 
Cts / 1 — M -2 = 0.494CT As r increases, the ordinate a of J_ decreases to a minimum at 
point I< 2 with r 2 = 0.0918, increases to a maximum at point K\ with t x = 0.609, and then 
decreases reaching the origin P at t = 2.25 with M = 6.2125. Points K\_ and A 2 divide 
J into three segments J n , n = 1,2,3. Similar to Figs. 2 and 4, curve J partitions the 
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domain of dependence V into V + and T>~ . Also shown are the partition of T> ± under the 
rules to be described in §3. 

2.2 The critical curve J in the £<r plane 

Now we study the properties of the image J in a meridian plane, the £<r plane with 
a > 0. With i related to £ by (1.12), (2.2) and (2.3) become the equations for the image 
J in the £cr plane, 

£ = C(t-r)/M-C J <7 = C{t — t)\J M 2 - 1/M (2.5) 

and 

df CM + RM' da r—p. , da . . 

— = \-CM = - — y/M 2 — l hence — = — tan a , (2.6) 

dr M 2 dr dt, 

Thus the tangent line to the curve J" is a backward Mach line, 

a - Rcosa = -tana{£ - [X(r) - X(t ) + Psina]} . (2.7) 

Note that the image of a semicircle R = C(t — r) in the £a plane is a semicircle of the 
same radius in the £cr plane centered at (— X(t) + A(r),0), 

K + X(t) - X(r)] 2 +a 2 = C\t - rf . (2.8) 

From (2.7) and (2.8), we see that 

{9} The curve 3 in the plane is the envelope of the family of Mach lines (2.7) and also 
of the family of circles (2.8) with parameter r e (t*,t*). 

Figure 3 shows the images of curves in the £a plane shown in Fig. 2. At a constant 
supersonic speed, the curve J becomes the backward Mach line TP enveloping the family 
of circles (2.8). Figure 5 shows the images of curves in Fig. 4 for a constant accelerating 
motion. The curve J begins at point T at r = 1 and M = 1 with a vertical tangent 
corresponding to the back facing Mach angle 7r — a = 7r/2. The curve approaches the 
origin P as r — + t along a radial line which is the backward Mach line. Tangent lines to 
J or the local backward Mach lines are shown in the insert. Figure 7 shows the images 
of curves in Fig. 6 for a constant rate acceleration from a low supersonic Mach number. 
The semicircle B with radius Ct and centered at (-A(f),0) is the image of the circular 
boundary B of f> when r = 0. The critical curve J is divided by the point 7\ , of local 
maximum and the point A '2 of local minimum into three branches. See also the insert. 
The first branch J\ starts at the origin P along the radial line or the the backward Mach 
line at t = t. The third branch J 3 ends at point T on the semicircle B with a common 
tangent; the backward Mach line at r = 0 with M = 1.15. 

Using (2.6), we obtain the curvature of J, 

Mj _ M'{t) (29) 

* M 2 a'(r)M 2 VAP - 1 
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The sign of k depends on the sign of M' j a' . From {6}, we know that at an extremum 
point on J , o’ = 0 , M' > 0 and (2.9) says that the radius of curvature 1 /k at the image 
point on J vanishes. If the extremum point on J is either a local maximum or minimum, 
say point K n at r„, then o' and 1 /k change sign along J as r passes r n . The point I\ n , 
which is the image of a local maximum or minimum point K n on jf , is a cusp on J formed 
by two adjacent branches J n and J n +i- 

Now we study the curve J or rather the two adjacent branches in the neighborhood 
of the cusp K n . From (2.3) and (2.9), we see that on these two branches, o' , and 1 /k 
are of opposite sign and they vanish at K n . With M' > 0, the slope of J is increasing 
continuously as r increases across r„ from branch J n to J n +\. We describe the curve J 
in the neighborhood of a cusp K n , 

The image of a local maximum (minimum) point K n on J is also a local maximum 
(minimum) point on J. The segments J n and J n +i on the left and right of point I< n 
are mapped onto the right and left branches J n and J n +i of the cusp K n . 

As r increases across r n , the slope of the tangent line or the backward Mach line (2.6) 
increases along branch J n +\ to the cusp K n and then continues along branch J n . If a 
attains a local maximum (minimum) at cusp K n , / reaches its minimum (maximum). 
The branch J n +i is concave downward (upward) and lies below (above) J n when 
the cusp K n is a local maximum (minimum) point. The two branches of opposite 
concavity form a simple cusp of the first species. 

From (2.3a) and (2.9), we know that if M' = 0, then o’ < 0 and k = 0. and that the 
curvature of a branch can change sign if M' changes sign. Also we see that M' has to be 
positive, when a' > 0. Therefore, 

{11a} A point on J7, where the Mach number is a local maximum or minimum, is a point of 
inflection of J . The curvature of J changes sign whenever the velocity changes from 
decelerating to accelerating or vice versa. 

{116} A branch joining a point of local minimum a to a local maximum, as r increases, has 
to be in an accelerating subregion and has to be concave downward. A decelerating 
interval, M' < 0, can occur only on a branch joining a local maximum to minimum 
and this branch is concave upward during an accelerating interval but downward in a 
decelerating interval. 

To see the additional partitions needed for the subdomains V' ± , so that the mapping 
(1.12) of each subdomain on to its image is one to one, we study the mapping of a horizontal 
segment S r S{ of constant <r c in the (a plane to its image S r Si in the £cr plane. If the segment 
S r Si lies either in V + or in V~, i. e., the segment does not intersect the critical curve J , the 
Jacobian J = (1.13), does not change sign and the mapping is one to one. We need 

to study only those segments crossing J , on which J = 0 or cos ip = \/M 2 — 1 /M > 0. 
Therefore, we consider the segments with £ > 0 that cut two adjacent branches of the 
critical curve J, the left branch J n and the right branch J n +i, at 5„ and 5 n+ i as shown 
in Fig. 8. In Fig. 8a (8b), the point K n is a local maximum (minimum) and o c is less 
(greater) than a n . The segment S r §i is divided into three segments 5 r 5„+i, 5 n +i5 n and 
S n Si. They are oriented as indicated by the arrow in the direction opposite to the <f axis as 
r increases, because = —CR/£ < 0. Since ±d£/d£ > 0 in V’ ± , the image of the segment 


{10a} 

{ 106 } 

{10c} 
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of S r Si in V ± is oriented in the direction of qF<£ axis as shown in Fig. 8. Note that the two 
image branches J n and J n +i form a cusp at K„ with branch J n lying to the right (left) 
of Jn+i- in Fig. 8a (8b). The first and third segments S r S n +i and S n Si lie above (below) 
J therefore, their images S r S n + 1 and S n Si are oriented in the same (opposite) direction 
from right to left (from left to right). The second segment S n +iS n lies below (above) J , 
therefore its image S n +iS n is oriented in the opposite (same) direction from right to left 
(from left to right). Because of the change of direction of the image segments across each 
branch, the second segment S n +iS n and the two branches J n and J n +\ lie on the same 
sheet while the first and third segments lie on two different sheets. Therefore, the first 
segment S r S n +i will be continued across the branch J n without a reversal of direction and 
the third one S n S[ across J n +i- 

Consider the inverse transformation of (1.12) of a point in the £ a plane without 
considering the three separate sheets. A point q in the middle segment S n +iS n has at 
least three image points q r ,q and qi which lie in the segments S r S n +i, S n +iS n and S n Si 
respectively. As point q approaches an end point say S n , the two image points qi and 
q approach each other and coincide at S n as a double root. The third image point q r 
approaches a point q* to the right of 5„+ 1 - When q crosses over to the right of 5„, there 
is only one image point qi which moves to the right of the point q*. We say that 

{12a} A horizontal segment S r Si intersecting J twice appoints S n and 5 n +i is divided by 
J into three nonoverlapping horizontal segments S r S n +i, S n +iS n and S n $i oriented 
in the same direction. They are mapped respectively onto three segments S T S n + 1 , 
Sn+iSn and S„Si with folding or change of orientation at the points S n and S n +i on 
the branches J n and J n +\ respectively. 

{126} The segments S r S n +h and S n Si overlap each other over the segment S n +iS n - A point 
S on Sn+iSn lies on all three segments; therefore, the multiplicity of point S is three if 
the inverse mapping is restricted to the segment S r S[. If point 5 lies on either S r S n + i 
or S n Si , the multiplicity of 5 is one. 

To end this section, we recall that a critical point K e on J which is neither a local 
minimum nor maximum is excluded in the above discussion because the image of 7i e is not 
a cusp and behaves the same as a point on J with o' ^ 0. If a parameter A (or parameters) 
defining U(t) can be varied so that two adjacent roots r n and r „+ 1 of o’ = 0 for A < A e 
become a double root r e at A = A e and disappear for A > A e , then the extremum point at 
the double root r e is neither a maximum nor a minimum. There is no cusp formation in 
the mapping of the neighborhood of point K e at r e when the parameter A > A e . A small 
decrease in A from A e results in two cusps appearing like a butterfly [4] in the mapping of 
the neighborhood of K e . There is an increase in multiplicity of 2. The sensitivity of the 
mapping on the parameter has been known in many other physical problems [5] and will 
not be addressed here. 

As a complement to the above quantitative analysis, we present in Appendix B the 
qualitative analysis of the inverse transformation of (1-12) and derive the formula for the 
retarded time for a point in the neighborhood of a branch near or away from a cusp. In 
§3, we use the above statements {1} to {12} to formulate the rules for the partition of the 
domain of dependence f>. 
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3. The Rules for the Partition of the Domain of Dependence 

Guided by the statements {1} to {12} in §2, we formulate the general rules for the 
partition of of the domain of dependence V for a medium moving at an unsteady velocity 
U(r)i. Since the multiplicity of V is 1, we seek to partition V into I nonoverlapping 
subdomains V such that the mapping from V to Vi in the £a plane by (1.12) is one to 
one. Hence the multiplicity of V t is 1, i. e., for each point in V, there is only one retarded 
time. Also we define the minimum number of partitions I. 

Since different supersonic intervals in [0, t] are separated by subsonic intervals, we can 
study the partition of the domain of dependence for each supersonic region one by one. 
Therefore, it suffices to study the case where there is only one supersonic interval 
The Jacobian J of transformation (1.12) changes sign in the supersonic region across the 
critical curve J on which J = 0. The curve ends at r = t* on the £ axis. It begins at 
r = f* from the £ axis if t* > 0 or from the semi-circular boundary B if = 0. Therefore, 
the critical curve J partitions V into two nonoverlapping regions V ± where ±J > 0. The 
region V~ is bounded below by the £ axis and above by the the critical curve J if t* > 0 
and in addition partially by B if t* = 0. From (2.1) we see that J and hence V" - lie inside 
the circular sector 0 < V < Vm where t/> m = arccos[l/Af moI ] and M mai is the maximum 
of M(r). 

From {12a} and {126}, we see that the mapping from V* to their images V ± may 
not be one to one. Further partitions of V ± are needed to separate the region in V + 
above ( V~ below) the curve J to the left hand side of a peak (valley), a local maximum 
(minimum) point K n at instant r n where n is odd (even), from the right hand side of the 
peak (valley). There is not a unique way to do the separation. We choose to partition V + 
(V -- ) by the large (small) circular arc (C~) passing through the point K n with radius 
C(t — r n ) at the retarded time r n for n odd (even). 

Before we state the rules for the partition of the domain of dependence, f>, for a given 
velocity function U(t), t 6 [0, f] with a supersonic interval, we carry out the following two 
preparatory steps. 

Step 1. Graph the boundary of the domain of dependence V in the £cr meridian plane, 
bounded by a semi-circle B with radius Ct centered at the origin, and then the critical 
curve J defined by (2.2). Locate the roots, of a' = 0 in (2.3a) while a" 7^ 0 and arrange 
the roots, r n , n = 1, 2 • • • , N in decreasing order of t, i. e., r„ > r n+ a . These N points 
divide J into N + 1 segments Jj , j = 1, • • • , N + 1. Segment Jj joins points Kj-i and 
Kj. Here Kq and AT/v+2 denote the end points of J at t = t* and r*. 

In this step we made use of {6} and {7} and note that these roots should be in accelerating 
supersonic interval(s). The first segment J\ joins the end point at r = /' on the £ axis 
to the local maximum point K\. When r* > 0, the last segment 3n+ 1 joins jKiint I\ N to 
the end point A'/v+i on the £ axis. Hence Kn is a local maximum and .V has to be odd. 
If r* = 0, the end point A'jv+2 lies on the semicircle B of r = 0. Kn is a local maximum 
(minimum) if N is odd (even) while N is odd (even) if er'(0) is positive (negative). Let 
L + and L~ denote the number of maximum and minimum points on J . then we have 
L + + L~ = N and 

L + = L~ + 1 = (N + l)/2 when N is odd, (3-1) 
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L + — L — N/2 when N is even. 


(3.2) 


Step 2. Graph the images of the boundary curve B and the critical curve J in the £a 
plane. The image B is a semi-circle of radius Ct centered at (-X(t), 0). The image J 
forms cusps at the N image points K n . 

These N cusps divide J into N + 1 branches, . the images of the segments J n . The 
domain od dependence V is the maximum area bounded above by either B or J and below 
by the £ axis. 

Now we state the rules of partition: 

Rule 1. Partition T> by the critical curve J into where ± J > 0 

The domain V~ is bounded below by the £ axis and above by J if the velocity begins 
subsonic M(0) < 1 and becomes sonic at r = t*. If the velocity is supersonic at r* = 0 
the end point T or Ky+\ lies on the semicircle B with polar angle V’iv+i = arccos 1/M(0) 
and then V~ is bounded above in addition by the small arc C^ +1 of the semicircle B 
with 0 < ip < *Pn+ i- Note that the large arc C^ +1 of B , with V’n+i < < ft, is a 

boundary of V + . A semicircle of constant r in the supersonic region is divided by the 
critical curve into two arcs C ± in V ± . The arc C~ is the smaller one since its central angle 
is arccos 1/M < 7 t/ 2. The subscript n will be added to C* when r = r„, n — 1, • • • ,1V 4- 1, 
with Tw+i equal to t*. 

Rule 2a. Partition V + into 1 + L + nonoverlapping domains v r +. 1 ,; = i,2,-.-,i + L+, 
by the L + large circular arcs, Cf l _ 1 , with constant retarded time T 21-1 passing through the 
local maximum point K 21-1 on J. 

Rule 2b. Partition V~ into 1 + L~ nonoverlapping domains V 2 J, l = l, 2, • ■ • , 1 + L~, by 
the L~ large circular arcs , C 2 {, with constant retarded time r 21 passing through the local 
maximum point K 21 on J . 

The subdomains and circular arcs in t> ± are numbered so that they can be grouped 
together as N + 2 subdomains V* and N + 1 circular arcs Cf at r,- using the plus or minus 
superscript for i odd or even. Thus the domain of independence V is partitioned into I 
nonoverlapping subdomains V±, where 

J = (1 + L + ) + (1 + L~) = 2 + N (3.3a) 


is the minimum number. 

A subdomain is bounded by the £ axis, two circular arcs Cf_ 2 , Cf and two 
segments of the critical curve and . Here, the segment J 0 reduces to the end 

point Kq of J and the arc Cf is empty when i < 0. 

Rule 3. Find the images of the I nonoverlapping subdomains in the £<7 meridiem 
plane one by one and denote the images by using the plus or the minus superscript for 
i odd or even. The multiplicity of a is one. 

Since the mapping (1.12) is independent of the radial and polar coordinates, cr and 6 , we 
can now define the volumes of integration of the integral in (1.16) in the £,77, ( space. 
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Rule 4. Rotate each one of the N + 2 subdomains Vj- in the meridian plane about the 
£ axis to generate a volume of revolution in £,??,£ space and denote the volume by the 
s am e symbol V* . The volume integral in (1.16) for the acoustic potential in the x space 
is expressed as the sum of the N + 2 integrals over the volumes V*. The multiplicity of 
each volume is one 

In case that there is more than one supersonic interval in [0,t], say s of them, we apply 
the above rules to these regions one by one. We locate the local maximum and minimum 
on the critical curves in the regions, draw the corresponding partition circular arcs in the 
£cr plane and find the images in the £cx plane. The total number I of nonoverlapping 
domains in the £cr plane or the number of domains of integration in the £a plane is 

I = 1 + 5 + N , (3.36) 


where N denotes the total number of maxima and minima on the s critical curves. When 
s = 1, we recover (3.3a). 

A formula or a numerical scheme is needed for the inverse solution of (1.12) to deter- 
mine the retarded time r for a point S in a V) . Consider a point So^o^o) in a domain 
V* defined under Rule 2 a or b. Let points Si and S r denote the left and right end points 
of the horizontal segment in the domain passing through point So- Each end point will 
either be on a circular arc of constant r or a branch of the critical curve on which o' ^ 0. 
Thus we know the values of r and £ and £ for both end points. Let us use the subscripts 
l and r to denote the values associated with points Si and S r . If point So coincides with 
one of the end points, the inverse solution of (1.12) is found. If So is not an end point, 
J = d£/d£ ^ 0 in the segment, hence there is a unique inverse solution £o- Now we 
can determine the numerical solution £o of (1.12) by Newton’s method using the linear 
interpolation formula from the end points to get the first estimate for £o . If the point So 
is close to a cusp of the region or to an end point on a branch of the critical curve, we 
sba.ll compute the first estimate of t 0 by the approximate formula for the roots of (1.12) 
presented in Appendix C choosing the root lying in between £j and £ r . Thus we have a 
procedure for the determination of the retarded time To and provide a quantitative proof 
that the multiplicity of point So in a domain V* is one. 

Now we are ready to express the integral in (1.16b) for the velocity potential at 
P{x, y, z) as the sum of integrals over the volumes of revolution V,, under Rule 4, 


W,z,y,z)=--L£j [j[ 


de g(T, X + £, y + o cos 0, z+_o sin 0) , ^ ^ 


|Ri — M,-£i| 


In (3.4), Vi denotes a subdomain in the meridian plane under Rule 3. For each Vi, r* and 
hence Ri , M, and £j are functions of £, o and t independent of 6. This formula does not 
need the knowledge of multiplicity. 

To demonstrate the meaning of multiplicity and the application of (3.4), we consider 
the source distribution to be a point source of unit strength created at point S(x', y',z ') for 
r > 0. If 5 lies in V^ , • • • , Vi m , with ij € [1, 1] and multiplicity m < I, then the intensity 
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received at P is 


(3.5) 


i(t,x,y,z) - - J^l/{4ir|fij ( - A'/,,?’,, I) • 

j~ 1 

Note that the domain of dependence V of point P(x , y, z) and the associated subdomains 
Vj are f-dependent. They were empty for t < 0, grow in size and change in shape as time 
t increases. Consequently, the multiplicity of point S is also f-dependent. The multiplicity 
changes by an integer and results in the addition or removal of term(s) on the right hand 
side of (3.5) and a discontinuity of the intensity at P. 

Now we formulate the rules for the change of multiplicity of point S, or point S(£, cr) in 
a meridian plane, when it enters or leaves a subdomain . Besides the axis of symmetry, 
the £ axis, a subdomain is bounded by a segment of the semi-circle B, branch(es) of the 
critical curve J and/or circular arc(s) Cf 1 , i < N , partitioning V ± . These three types of 
boundaries, which are moving, are called wave fronts. 

The multiplicity m(t) of point S begins with zero at t = 0. If point S remains outside 
of V for all f, we have m(t) = 0. The multiplicity changes when point S crosses a wave 
front in entering the domain of dependence V or rather a subdomain. The rules for crossing 
the three types of boundaries are different and hence are stated separately. 

Rule B. When S crosses B to enter (leave) the subdomain V + or V~ ofV, the multiplicity 
m increases (decreases ) by 1. 

Rule J. When S crosses a branch of J with cr'(r) > 0 from the left to the right hand 
side, the multiplicity m increases by 2, and with a' < 0, m decreases by 2. When the 
direction of crossing is reversed the change of m is reversed. 

Rule C. When S crosses a circular arc, Cf, i < N, partitioning V ± , there is no change of 
multiplicity, because the Jacobian J does not change sign and S is entering one subdomain 
while leaving another. 

These three rules can be applied to determine the multiplicity of point 5(£i,<Ji) in V 
at an instant f, by moving from a point S'o(£o ? ^o) with known mo, along a path to S\ and 
counting the changes of multiplicity along the path. For example, we can pick a point So 
outside of V with <r 0 = <J\ and hence m 0 = 0 and the path to S\ is horizontal. 

In the next section we shall demonstrate the applications of the rules of partition and 
counting of multiplicity for various unsteady velocities U(t). 

4. Examples 

For the cases where the velocity U(t)i of the medium relative to the coordinate system 
x remains subsonic [2], [3], the Jacobian of the transformation of V to V remains positive, 
and the multiplicity of the domain of dependence V in x is equal to one. These cases are 
examples of the class 1=1, where I is the minimum number of nonoverlapping domains 
in V. In the following, we present examples having maximum multiplicity greater than or 
equal to 2. These examples must have a critical curve, where J — 0, and hence a supersonic 
interval, in (0,f). See {2}. 

It was stated in {1} and {2}, that the domain of dependence V in the coordinate 
system x, with the medium at rest, is independent of the velocity of the medium relative to 
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the system x, the critical curve J is dependent on the velocity. The parametric equations 
for J in a meridian £cr plane are given by (2.2). The first step is to locate the local 
maximum and minimum points on the curve J , i. e., to locate the roots of ct'(t) = 0, 
(2.3a), with a" ^ 0, which are the roots of 


F(t) = (f - t)M' - M(M 2 - 1) = 0 , r € (0, t) , (4.1) 

with 

F\t) = (t- t)M" - 3 M 2 M' ± 0 . 

Here we set C = 1 and hence M = U(t). To demonstrate the rules of partition and 
counting in §3, we present first in §4.1 to §4.5 examples with only one supersonic interval 
in [0, <] i. e., 5 = 1, and then in §4.6 and §4.7 examples with more than one supersonic 
interval, s > 2. In the last two subsections there are also examples showing that the 
maximum multiplicity in V can be less than the number of partitions or the number 
of domains of integrations in (3.4), i. e., m < I. With s > 1 and N > 0, we have 
7=l + s + iV>2. We begin in §4.1 with 1 = 2. 

4.1. Class 1 = 2, the critical curve J has no local maximum. From {7}, we see 
that the curve does not have a local maximum or minimum. Thus (4.1) has no root, i. e., 
N = 0. The simplest example is that of constant supersonic speed [3], V = M' = 0 and 
M > 1. In the £cr plane, the critical curve J is a radial line PT which partitions the 
domain T> into V+ and V~ under Rule 1. Since J has neither a local maximum nor 
minimum, further partition of V’ ± is not needed. In Fig. 2, the coordinates a and £ are 
scaled by Ct and the concentric circles are circles of constant r/t with radius 1 — (r/<). 
The radius of the boundary B is 1. Thus Fig. 2 remains the same for allt > 0. 

The mapping of Fig. 2 onto the £cr plane scaled by Ct will also be invariant for all 
t > 0, as shown in Fig. 3. The superposition of the images V + and V~ yields the domain 
of dependence T > , i. e., V + U V~. The domain T> is bounded above by the critical curve 
PT, a backward Mach line, and by the large circular axe QT of the semicircle B at r = 0 
and below by the £ axis. Also shown are the semicircles at retarded time r 6 (0, t). These 
semicircles are divided by their envelope PT into smaller and larger arcs. The family of 
the larger (smaller) arcs covers V + (V - ). Thus V = V + and the multiplicity of T>(1 V~ is 
two and that of T> ^ V~ is one. Detailed discussion of this partition can be found in [2]. 

It is well known [3] that the domain of dependence of a point P(x, y, z ) with the 
medium moving at a constant supersonic speed is the backward Mach cone from P. Since 
the source distribution has been turned on all the time, the retarded time is in (— oo,<]. 
The multiplicity for a point in the cone is always 2. In the current initial value problem 
the source distribution is turned on at r = 0, therefore, the retarded time has to be in 
[0, t] and the domain of dependence shown in Fig. 3 is a finite part of the backward cone 
with multiplicity 1 or 2, as defined in the preceding paragraph. 

We use this simple example to show the meaning of multiplicity to an observer at P. 
We consider a point source at S of unit strength initiated at r = 0. Since both points 
P and S are fixed in x, S is a fixed point in the relative coordinates £,y,C and in the 
£, a plane. But in Fig. 3, the position of point 5 is moving along the radial line from 
point 5(1) towards the origin P (shown as the dotted line) with the instantaneous radial 
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distance |5(f)P| equal to |5(l)P|/f where 5(1) is the position at t = 1. The point 5(1) 
shown in Fig. 3 is lying inside the backward Mach cone but outside V. That is, at t = 1 the 
multiplicity of point 5 is zero, or the signal from 5 has not reached P. At an instant t a > 1, 
the radial line intersects the larger circular arc QT of B and enters the region V V~ with 
multiplicity 1. At a later instant tb > t a , the radial line intersects the smaller circular arc 
of B and enters the region V~ with multiplicity 2. The observer receives the signal from 
5(fi) with intensity l/(47ri?*) when tb > t\ > t a and with the intensity doubled 1/(2ttR*) 
when t > tb- Here R* denotes the pseudo-distance (1.17). 

Additional examples for this class of I = 2 will be found as “degenerated” cases of 
the classes of I > 2 in the following subsections. 

4.2. Class 1 = 3, the critical curve has one local maximum. The simplest example 
for this class is the case of constant acceleration from zero to supersonic speed. By choosing 
the time scale, we can write M = r and then the the supersonic interval is from i* = 1 to 
t* —f > 1. Equation (4.1) has one root, r x = t 1 / 3 . The curve J begins from the £ axis at 
r = 1 with M* = 1 rises to the maximum K\ at r x = t 1 / 3 and descends to the origin at 
r = t. Thus this case belongs to the class of N = 1 and I = N + 2 = 3. Note that we can 
no longer get a time invariant critical curve and partition of the domain of dependence by 
using Ct as the length scale and t as the scale for r as we did for constant supersonic speed 
in §4.1. We note that when t < 1, M < 1 for r € [0, <], J does not exist. When t > 1, 
M > 1, J exists in the interval [1, t] but J and the interval do not grow proportionately 
with t. We have to construct the critical curve and the partition for a given t. 

Figure 4 shows the critical curve J for t = y/E. Besides applying Rule 1 to partition 

V into V" + and V~ , we have to apply Rule 2a to partition V + by the circular arc C+ of 
constant r x passing through point K\ into V"* and V+ . Thus V is partitioned into 7 = 3 
nonoverlapping domains, Vj + , F 3 + and V~ . 

The images of the nonoverlapping domains in Fig. 4 are the three domains of integra- 
tion in £ and a in (3.4). The superposition of these three domains for t = y/E onto the £<r 
plane is shown in Fig. 5. The union of these three images gives the domain of dependence 

V which is bounded above by the semicircle B and by branch Ji and below by the £ axis. 
Besides the £-axis, V~ is bounded by the two branches J\ and J 2 , V 1 + by the circular arc 
C+ and the branch J\ and V 3 + by the semicircle B, the arc Cf and the branch Note 
that V~ € Vj + , fl Vj + € V~ and V~ D F 3 + ^ 0. The multiplicity of V~ Pi V 3 + is 3, 
that of y- i y 3 + is 2 and that of T> £ V is 1. More detailed descriptions of these three 
domains of integrations are given in [2]. 

We use this example to explain that the domain of dependence V, the semicircular 
disc of radius t in the meridian plane, is the domain for all source fields q(r,x,y,z ) in 
the interval (0,t). But for a given source field, there is an effective domain of dependence 
which is a subdomain t>. For example, if the source field is created at a time b > 0, i. e., 
q = 0 for r < b, then only the part of V where the retarded time r lies in the interval 
[b,t] will contribute to the acoustic field at point P at time t. This effective domain of 
dependence V e is the subdomain of V below the semicircle B e of r = 6_or radius R = t — b. 
Thus only the portion of the critical curve J and the subdomains lying inside this 
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semicircular disc t> e will be mapped to the (a plane. f For the current example, M(r) = r, 

(4.1) becomes 

F(t) = t — t 3 = 0 , T € (M) , ( 4 - 2 ) 

with F'(t) = -3 t 2 < 0. Equation (4.2) has one root r x = t 1 ? 3 provided that b < t 1 / 3 . 

For b < 1, the supersonic interval (l,t) is in (b,t). Hence the critical curve J and the 
large circular arc C* remain inside V e . The subdomains V~ and Vj + shown in Fig. 4 remain 
unchanged. Only the subdomain U 3 + is reduced in size with the semicircular boundary B 
replaced by a smaller one B e . Consequently, the images J, V~ and V+ in the £<r plane 
remain the same as those shown in Fig. 5. The image F 3 + is reduced in size with the 
boundary B replaced by B e . 

If n > b > 1, part of the segment J 2 and the subdomain V 3 + in addition to that of 
V~ in Fig. 4 are outside of T> e or cut off by the semicircle B e and their images will be cut 
off in Fig. 5. Only the image Vj + remains unchanged. 

When b = T\ , the semicircle B e passes through point K \ , coincides with the semicircle 
C\ at T\ and completely cuts off the segment j 2 and the subdomain U 3 + in Fig. 4. Thus 
the critical curve J e = J\ does not have a local maximum. The domain of dependence V e 
in the £cr plane is equal to V+ in Fig. 5. The image V~ is bounded by the £-axis, J 2 and 
the small circular arc Cf which together with the arc C~l form the semicircle B e . 

When t > b > T\, the semicircle B e not only cuts off the entire branch J 2 , the circular 
arc Ci and the domain U 3 + but also parts of the branch J\ and the domains V~ and V+ . 
The critical curve which is the remaining part of J\, decreases monotonically to the 
£ axis as r increase from b to t, i. e., (4.2) has no root. Thus our example degenerates to 
the class of I = 2 when t > b> T\. 

For a fixed b > 1, (4.2) defines a critical time t c = b 3 for which F(b ) = 0, t c = T\ and 

(4.2) has no root. As t increases over # c , I jumps from 2 to 3. That is, an observer at point 
P will receive signals from point 5, fixed in the coordinate system x, twice if 5 € V~ when 
t < t c with 1 = 2 but three times if 5 happens to be in V~ U F 3 + when t > t c with 1 = 3. 

4.3. The class 1 = 4. The critical curve J has one local minimum and one 
maximum. We consider the case of constant rate of acceleration from an initial Mach 
number M 0 , i. e., M(r) = M 0 + r 2 , with F(r) = 2 r(t — r) — M(M 2 — 1) and F" = 
—2 — 6M(M + 4r 2 ) < 0. 

If Mq < 1, the supersonic interval begins at r = > 0 with M* = 1, <j* = 0 and 

> 0. Using {8} and F" < 0, (4.1) has only one root and the curve J has only one 
maximum point. Thus the case for Mo < 1 belongs to the class 1 = 3. 

To create an example of the class I = 4, we consider Mo > 1 with f* = 0 and t* = t. 
With F( 0) < 0 and F(t) < 0, or a'(0) < 0 <r'(f) < 0 , (4.1) can have two roots or none. It 
is clear that M$ - 1 has to be sufficiently small so that F(r) can be nonnegative in (0, t), 
i. e., (4.1) has two roots. For a given t, there is an upper bound below which two distinct 
roots exist. When Mo is equal to the upper bound M c o, these two roots coincide in one 
double root r c for F(t c ) = F'(r c ) = 0. This and the equation M c = M c0 + r 2 , define r c , M c 


f This problem can be identified with the one with the source field created at f = 0 
while M = b by the time shift, f = t — b and t = t- b with velocity U(t) = U (r). 
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and r c o as functions of f.f For example when t = 2.25, we have M c = 1.382, t c = 0.3328 
and M c o — 1.275. 

If the initial Mach number Mo > M c o, (4.1) does not have two distinct roots and 
hence the critical curve J has no local maximum. Hence this case belongs to the class of 
1 = 2 in §4.1. 

For 1 < Mo < M c o, (4.1) has two distinct roots, the critical curve J has a local 
maximum and a local minimum as shown in Fig. 6 and the image J has two cusps forming 
a ’’butterfly” as shown in the insert in Fig.7. In those two figures, we have M(r) = 1.15+r 2 , 
r £ [0,2.25], i. e., Mo = 1.15 and t = 2.25. Since Mo < M c0 , (4.1) has two distinct roots 
n = 0.6093 and r 2 = 0.0918. We have N = 2 and I = 4. 

As shown in Fig. 6, the critical curve J begins with <r(0) = 1.1111 at point T on the 
semicircle B, decreases to a minimum at point A' 2 along the segment Jz, increases to the 
maximum at point K\ along J 2 and then decreases to the origin P along J\ . By Rule 1 
the critical curve J partitions the domain of dependence V into V ± . By Rule 2 and 3, 
these two domains, V ± , are partitioned respectively by the large circular arc Cf passing 
through K\ and the small circular arc C 2 passing through A 2 into V’ 1 + , Vy" and V 2 ~, VJ . 
Thus the domain of dependence V is partitioned into four nonoverlapping domains. 

Figure 7 shows the images of the critical curve, the semicircle when r = 0 and the 
two circular arcs in the £< 7 plane. Because the initial Mach number M(0) is above 1 and 
below the threshold value M c 0 , the circular arc C 2 is approaching the semicircle B from 
the right hand side while the branches J 2 and Jz and the two cusps K\ and K 2 can be 
seen only in the insert with 20 x enlargement. The domain VJ bounded by B,C 2 Jz and 
the £-axis is disappearing. In addition to the £-axis, V 2 ~ is bounded by J 2 , J\ and C 2 , V+ 
by Cj and J\, and V 3 + by B. Jz, Jz and C+ . The branch Jz begins at point T tangent to 
the circle B, goes downwards on the right hand side of B as the envelope of the circles of 
constant r £ [0.r 2 and ends at the cusp J\' 2 . The branches Jz and Jz and the semicircle B 
form a curvilinear triangle which is the intersection of the four domains of integration in 
the £cr plane. The intersection has the maximum multiplicity 4. 

If we increase M(0) to the threshold value M c 0 , the two cusps I\\ and J\ 2 and the 
branch Jz disappear. Branch J\ joins smoothly with Jz to form a smooth critical curve 
J with no local maximum or cusp. The partitions by the arcs C+ and C 2 are no longer 
needed. Thus we recover the class of I = 2. If we decrease M(0) to 1, the minimum point 
K 2 moves down to the £-axis and leftward to B while the branch Jz and the domain VJ 
disappear. Thus we recover the class of I = 3. 

Similar to the discussion in §4.2, for a given M(0) = 1.15, the critical time of observa- 
tion, when M c o = M(0), is t c = 1.074. When t < t c , an observer will receive signals from 
point S at most twice in the class of I = 2. When t exceeds t c , the observer can receive 
signals from the same point 5 four times in the class of / = 4. 


f From F'(t c ) = 0 we get r c 2 = (M 2 — 1)/(6M C ) and from F(t c ) = 0, t 2 = (M 2 — 
1)(3M 2 + 1 ) 2 /(6M c ). The second equation defines t > 0 for a M c > 1. The inverse 
transformation is also unique. Thus for a given t > 0, the second equation yields a unique 
M c > 1 and then the first one, r c > 0. They in turn define the corresponding initial Mach 
number M c 0 = (5 M 2 + 1 )/6M c > 1. 
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4.4. I — 5. The critical curve J has one local minimum and two local maximum. 

It follows from {7} and {8} that K 2 is the local minimum point and K\ and A' 3 are the 
two maximum points, with T\ < t 2 < t 3 and that > 0 and A/(f*) = 1 for t* > 0 or 

M(U) > 1 for U = 0. 

We consider the case that the acceleration is a quadratic function of t, 

M'(r) = c 0 + 2cir + 3c 2 t 2 , (4.3) 

and hence M is a cubic function, 

M(r) = [c 0 +C 1 T + c 2 t 2 ]t , (4.4) 

with co + ci +c 2 = 1. We have chosen the time scale such that M( 0) = 0 and M( 1) = 1. As 
explained in §4.2, a nonzero initial value M(0) can be converted to (4.4) by an appropriate 
time shift and change of coefficients Co and c\ . 

We consider the case that the speed U = M is accelerated to supersonic M a > 1 at 
r = a > 1, followed by an interval of deceleration to Mb at r — b and then acceleration to 
r = t. To show the accelerating and decelerating intervals explicitly, we rewrite (4.3) asf 

M 1 = c(r - a)(r — b ) with b > a > 1 and c = 6/[2 — 3(a + b) + 6ab] . (4.5) 

To insure that M > 1 during deceleration, we assume that M(b) > 1. 

Using (4.5), we find F( 1) > 0, F(a) < 0 and F'(r) < 0 in the interval (1, a), therefore, 
(4.1) has only one root in the interval. Since F(b) < 0, F(t) < 0 while F" (t ) < 0 in the 
interval (&,<), (4.1) has either two roots or none in the interval (b,t). Therefore, J will 
have a maxim um in the first acceleration interval, (1, a) and a minimum and maximum or 
none in the second acceleration interval, (b,t), in agreement with {8a}. Whether there is 
a pair or not in the second interval, i. e., whether I = 5 or 3, depends on the values of a, b 
and t. Let us consider the function M(r) specified, i. e., a and b specified while t varies, 
as we study the signal received at point P from a point S. It is clear that when 1 < t < b, 
(4.1) has only one root, and J has only a maximum. We have the class of / = 3. There 
is critical time t c > b when there is a double root Td of (4.1) in the interval ( b,t c ). By 
eliminating (t c - r d ) from F{r d ) = 0 and F'(T d ) = 0, we get an equation for r d , 

M"(r d )[M 2 (r d ) - 1] - M(T d )[M'(T d )} 2 = 0 (4.6) 

and then compute t c from (4.1). We have an example for the class of I = 5, when t > t c . 

Figure 9 shows the domain V and the critical curve J for a = 1.7, b = 3^05 and 
t = 5. Now t c = 3.666 < t, we have an example for the class of I = 5. The curve J which 
partitions V into V ± , has a maximum point Ah (0.880, 0.642) at Ti = 3.911, a minimum 
point A'2(l-941, 0.287) at r 2 = 3.077 > b and then a maximum point A” 3 (3.078, 1.641) at 
t 3 = 1.511. The larger circular arcs C+ of T\ and of r 3 partition the domain V + into 

f Here we assume that the quadratic function (4.4) has two roots in the interval (l,f). 
We can have examples of I = 5 with M' = [(r + a) 2 + 6 2 ]/[( 1 + a) 2 + b 2 ] > 0. See the 
report [6] 
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three, F 1 + , V 3 + and F 5 + . The smaller circular arc C 2 partitions the domain V into two, 
V-andV". 

The images of these five nonoverlapping subdomains of T> in the £<7 plane are defined 
one by one and the acoustic potential is given by (3.5) as the sum of integrals 

over the five images V 1 + , V 2 ~, F 3 + , VJ and F 5 + . These five image subdomains are super- 
imposed onto the £a plane as shown in Fig. 10a. Their union is the domain of dependence 
V. The five subdomains Vj~ can in turn be identified from Fig. 10a by tracing the images 
of the boundaries of tq . 

Fig. 10b shows the critical curve J enlarged five times in order to display the four 
branches and the curvilinear triangle bounded by the three branches J \ , Jz and Jz which is 
the intersection of the five subdomains. The critical curve and the part of the semicircular 
B in the enlargement are shown in heavier lines than those for the three circular arcs Cf , 
C 2 and because multiplicity does not change across these three arcs, according to Rule 
C. The numerals in bold face indicate the multiplicities in a region 7Z bounded by B, J 
and/or the £ axis. We demonstrate the counting of multiplicities by applying Rules B 
an d J instead of identifying 71 as the intersection of subdomains V* and counting the 
number of subdomains in the intersection. For example, we show by counting that in the 
curvilinear triangle the multiplicity is 5. 

Consider the change of multiplicities from left to right along a horizontal line C of 
constant <r c below the cusp K\ and above K 2 , i. e., <7 2 < <x c < oi. A point on C lying to 
the left of the semicircle B is outside of V and has multiplicity m = 0. When C crosses 
B from left to right and enters V 5 + , the multiplicity m changes to 1 according to Rule 
B. When C crosses the branch J 4 , m = 3 according to Rule J . When C crosses J 2 and 
enters the curvilinear triangle, m = 5. When C leaves the triangle and crosses J\ . Jz (or 
Jz, Ji), and then B, m = 3, 1 and then 0. If it crosses Jz , B and then J\, m = 3,2 and 
then 0. If the horizontal line is below I\ 2 , m changes from 0, 1,3,2 and then 0. If the line 
is between K\ and Kz, ra = 0, 1, 3, 1 and then 0. If the line is above Kz but below <7 = t, 
m = 0, 1 and then 0. 

It is clear that we can find examples which require an even greater number of par- 
titions, but we have given plenty of examples to demonstrate the rules of partitions and 
counting of multiplicities. We note that although the velocity U or Mach number M in 
the examples are polynomials of r, the rules are applicable when M is not a polynomial. 
This is demonstrated in the following example. 

4.5. The speed U(t) is periodic We consider JJ = M = uq — (c/a)cos(ar + b ) and 
t = 7r /a. Fig. 11 shows the domain of dependence T> and the critical curve J with a — 
1,6 = 0 and c = 0.8a and uo = Mq = 2 and t = tt. The speed is supersonic in the interval 
[0,<]. The critical curve begins on the circle B and ends at the origin. There b one local 
minimum K 2 and one maximum Ki at r equal to 0.28194 and 0.67991. With .V = 2 and 
I = 4, the partition of V into four subdomains is similar to that in Fig G Figure 12a 
shows the image of Fig. 11 in the £<7 plane. Figure 12b shows the 45 x enlargement of the 
region containing the cusps K\ and K 2 and the area bounded by the three branches, J \ , 
J 21 and Jz , only in which the multiplicity is the maximum, 4. Again note the similarity 
of Fig. 12 to Fig. 7. 

In the following subsections we shall present examples with more than one supersonic 
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interval in [0, <]. Let there be two intervals ( 42 , 43 ) and ( 44 ,^ 3 ) with 0 < 44 < 43 < 4 2 < 
l 1 < t. The speed is subsonic in ( 43 , 42 ) and sonic at the end points. From {7}, we see 
that there is at least one maximum on the critical curve in the first interval (t 2 , 43 ). In the 
second interval, the critical curve may not have an extremum, if = 0 and the speed is 
decelerating from supersonic speed M( 0) to sonic speed. Since N > 1 and I = l + s+N > 4 
for s > 2, we deal first with the class 1=4. 

4.6. Class I = 4, s = 2, two critical curves having only one local maximum. As 

a simple example, we consider a linearly accelerating motion in the interval [0, t], 

M'(r) = T — b, with &€(0,t). (4-7) 

The velocity function is 

M(t ) = Mo — br + t 2 /2 , with Mo > 1 and M(t ) > 1 . 4.8 

We choose M 0 , b and t such that the minimum Mb = M(b ) = M 0 — b 2 j 2 is subsonic i. e., 
1 > Mb > 0. The speed becomes sonic at 12 and * 3 , where 

<■2 = b + \/b 2 — 2c and 63 = b — y/b 2 — 2c , (4-9) 

with c = Mo — 1. Here we have two supersonic intervals { 12 , t) and [0, 43 ). Only the critical 
curve in the first interval has a maximum at point K\ at the retarded time T\ € ( 12 , t). 
As shown in Fig. 13, the domain of dependence V is partitioned into four nonoverlapping 
regions V’ 1 + , F 2 “, V 3 + and V 4 ~. The region V 2 ~ is bounded by the £ axis and the two 
segments of the critical curve, as r decreases from t to T\ and from ri to l 2 . The 

region V 4 ~ is bounded by the <f axis, the small circular arc B~ and the critical curve J 2 
from 13 to 0. The region, where J > 0 , is partitioned by the large circular arc Cf of 
constant retarded time t\ into V’ 1 + and V” 3 + . The segments on the £ axis bounding V 2 ~ and 
F-, where J < 0 , are £ 1^2 and ^ 3 ( 4 . With £j = C(t — £>), we have (j > Here we 

have t\ = t , = 0 and L 4 = 0 , £4 = Ct. 

If the images of V 2 ~ and V 4 ~ do not intersect, i. e., V 2 D V 4 — 0, the maximum 
multiplcity of D can only be 3 < I = 4. This requires that the image segments £ 1^2 and 
£ 3^4 do not overlap. Noting the sign of J, we have (2 > 6 and £4 > £3 while £2 < £ 3 - For 
the two segments to be nonoverlapping, it is necessary that £1 < £ 4 , i. e., 

r l 1 [ ll M(r)dr 

^4 — £1 > / M(r)dr , or < M > = — < 1 . (4-10) 

Ji 4 l i ~ l 4 

This says that the average Mach number < M > or speed < U > in the interval [ 44 , 43 ], 
has to be subsonic. This statement is valid for any M(r) having two supersonic intervals 
and an intermediate subsonic interval in [ 44 , 43 ]. 

In Fig. 13, we have b = 1.0, Mo = 1.18, X(0) = 0 and t = 2.4. The speed is sonic at 
4 2 = 1-8 43 = 0 . 2 , Mb = 0.68 and < M >= 0.94 < 1 . The images of the four subdomains 
in Fig. 13 are shown in Fig. 14a. A 3x enlargement of the region containing the critical 
curves is shown in Fig. 14b where the bold faced numerals denote the multiplicity. In this 
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case the subdomains V 2 and V 4 do not overlap and the maximum multiplicity is equal 
to 3 < I = 4. 

Figures 15 and 16a and b show a counter example of Figs. 13 and 14a and b with the 
same t = 2.4, and smaller 6 = 0.8. The speed is sonic at l % — 1.33 and 13 = 0.27, with 
Mb = 0.70. There is a longer second supersonic interval and the average Mach number 
< M >= 1.18 is supersonic. As shown in Fig. 16, the subdomains V 2 ~ and V 4 do overlap 
and the maximum multiplicity is equal to I = 4. 

4.7.More than two critical curves, each having only one local maximum We 

consider an example with periodic velocity variation about a mean subsonic speed. Then 
there are three subsonic and three supersonic intervals. We choose A'(0) = 0, t = 5.0 and 
M'(r) = 0.57T sin(7rr) and M(r) = 0.75 — 0.5cos(7rr). 

As shown in Fig. 17, the domain of dependence T> is partitioned into seven nonover- 
lapping regions by the three nonintersecting critical curves i = 1,2,3 and the large 
circular arcs C l passing through the maxima A' 1 . Therefore J = 3 + 3-|-l = 7. The images 
of the seven subdomains in Fig. 17 are shown in Fig. 18a. A 2x enlargement of the region 
containing the three critical curves is shown in Fig. 18b where the bold faced numerals 
denote the multiplicity. In this case the maximum multiplicity is equal to 3 < / = 7. 


5. Conclusion 

In the interaction of an acoustic field with a moving airframe we encounter a canonical 
initial value problem for an acoustic field induced by an unsteady source distribution, q(t, x) 
with q = 0 for t < 0, in a medium moving with a uniform unsteady velocity U(t)i in the 
coordinate system x fixed on the airframe. Signals issued from a source point S in the 
domain of dependence V of an observation point P at time t may arrive at point P more 
than once corresponding to different retarded times, r in the interval [0, <]. The number of 
arrivals is called the multiplicity of the point S. The multiplicity equals 1 if the velocity U 
remains subsonic and can be greater when U becomes supersonic. For an unsteady uniform 
flow U(t)i, the acoustic potential $(t,x,y,z) is governed by the convective wave equation 
with variable coefficients, the velocity U(t) and the accelerati on U(t). The solution $(f,x) 
is given indirectly via the solution of the corresponding acoustic field «§ in the coordinate 
system x with the medium at rest. The solution (1.16) induced by the source distribution 
q(t,x) is 


Q(t,z,y,z) = J J J d£dr)d( 


q{r,x + Z,y + ri,z + () 


\R-Mt\ 

q(Tj,x + t,y + Ti',z + C) 


(5.1a) 

(5.16) 


To evaluate the integral over the domain of dependence V in (5.1a) we need to relate r, 
R and £ to t and the coordinates of the source point S(x + £, y + 77, z + £). The implicit 
relationships are: 


e = e + / u(t')dt' 


and R 2 = C(t — r) 2 = + y 2 + £ 2 . 


(5.2) 
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They yield multiple retarded times in [0,<] and the corresponding R t and /*. We need 
to find the roots of (5.2) and identify the domain Vi in the /, 77, / space for each r,-. Then 
<J> is written as the sum of integrals over the smallest number I of domains Vi in (5.1b). 
Here I denotes the maximum multiplicity of V. It is the purpose of this paper to formulate 
the rules identifying these dom ains Vi for a medium moving at an unsteady speed, U(t). 
First we study the Jacobian of the transformation x to x in §2 and then use these studies 
to formulate in §3 the rules for the partition of the domain of dependence, V in the x 
space, into nonoverlapping subdomains V) for an unsteady velocity, U(r) ,r £ [0, <]. We 
show that the partition has the required property that the value of r and hence R and / 
defined by (5.2) for a point (/,»?,/) in an image Vj is uniquely defined. Thus the domain 
of integration in (5.1b) is the sum of those images Vj' s. Since those images may overlap 
in the /,t/,/ space, signals from a point S in the space lying in m of those images will 
reach point P m- times. Hence m is th e multiplicity of point S. The rules for counting 
the multiplicity of a point in D are also stated in §3. The analyses in §2 and §3 justify 
the assertions made in [2] for the case of a constantly accelerating speed. Additional 
examples with unsteady velocities, accelerating and/or decelerating at different ranges of 
Mach n um bers are presented in §4 to demonstrate the applicatio n of the rules of partition 
of V and to show the singularities of J in the /, f?, ( space. 

Appendix A. Domain of dependence and the region of influence in a constant 
supersonic flow We study the domain of dependence V of an observation point P at an 
instant t > 0 for an acoustic field created in the interval [0, t] and the region of influence of 
a source point S during the interval [t, /] where t denotes the retarded time. Consider the 
meridian £<7 plane in the coordinate system /, 77, / containing points P and S. The /-axis 
is in the direction of the flow, i. Point P is the origin and V is bounded by the / axis, 
the backward Mach line from P and the large circular arc of radius Rq — Ct with r = 0 
and centered at G on the /-axis with / = -MR 0 . Note that the region of influence of a 
source point S for the interval [0, t] is the mirror image of the the domain of dependence 
of P with respect to the a axis and point S at the origin. Fig. 19 a and b show the two 
domains of integration V + and V~ for M = V 2 in a meridian plane. 

The backward Mach line from P is the envelope of the semi-circles C of constant 
t £ [0,f] with radii R = C(t - r). The Mach line divides a semi-circle into two arcs, 
the larger arc C + is convex while the smaller arc C~ is concave from the Mach cone. The 
domain V + coincides with V and is covered by the family of convex circular arcs C + . The 
domain V~ is covered by the family of concave circular arcs, C~ with V + V~ equal to 
the semi-circular disc of radius Ct centered at G. Consider a source point S in V~ £ V + . 
Point S lies on the intersection of two circular arcs in V ± , the convex circular arc C + 
of radius R + centered at G + with / = -MR + in V + and the concave circular arc C~ of 
radius R~ centered at G~ with / = -MR~ in V~. Note that these two arcs are associated 
with different semi-circles or different retarded times r*, with r + > r”. In Fig. 19a and 
19b, we choose r + = 0.75f and r~ = 0.25 1. The signal created at the retarded time t~ at 
point 5 in Fig. 19b propagates at time t onto a sphere of radius R~ while its center travels 
parallel to % to S~ with \SS~\ = MR~ = The region of influence of point S with 

retarded time t~ (in the meridian plane) is bounded by the two forward Mach lines from 
S and the circular arc of radius R~ centered at S~. Since points 5, S~, P and G~ are the 
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vertices of a parallelgram, point P lies on the small circular arc. The signal initiated at 
time t~ from point S reaches point P at time t. Likewise we see in Fig. 19a how the signal 
from S at the retarded time r + reaches P at time t. Similarly, we can see for the case of 
an unsteady moving medium the connection between the region of influence of point 5 in 
the interval [r, t] and the propagation of the signal from S initiated at r reaching point P 
at time t. 


Appendix B. Transformation of a Neighborhood of J in the £cr plane onto the 
£<r plane and the Inverse Transformation 

We study the transformation of a point S(£,cr c ) in the neighborhood of the critical 
curve J onto a point S(£,<t c ) near J and its inverse transformation. We add the subscript 
c to a to shown that a remains unchanged in the transformation. 

First we consider the point S near a branch of J but not too close to the extremum 
point K\. Using the same ordinate cr c for the neighboring point 5 c (£ c ,cr c ) on J\ , we use 
(1.9) and the equation for a in (2.2) or (2.5) to find r c or R c first and then £ c from (2.2), 
and £ c from (2.5), for the corresponding point 5 C (£ C , cr c ) on J . Here the subscript c denotes 
the value at r c . 

Note that we are considering, 

|£ — £ c | <C 1 and hence |r — r c | <C 1 while cr' c ^ 0 (J5.1) 

Here r denotes the retarded time corresponding to point S by (1.9), 

t = t - R/C and R = yj? + a\ . (B. 2) 

Using (B.l) we simplify the transformation (1-12) from £ to £. Using (1.12) and expressing 
M by the Taylor series in (r — t c ), we have 


Z-i, = i-L+J'cM( f) df = (-Z c +CM c (T-T c )+jM' c (T-T c f+0((T-T c f) , (B. 3) 
Now we express C[t — r c ) in power series of 8 C = (£ — £ c ), 


C(t - Tc) = ~{R - Rc) = -{d-^R) c 6c - -( d 2 u R) c 8 2 c + 0(S 3 C ) 


— —8c 

M c 




(BA) 


Using (2.3a) and (B.4) in (B.3), we get 


£-£c 


< + 0(8l) 


KRc - CAUM I - 1 ) f2 , 

2 CM?R c 

sfK^<n z 


2CRr 


'it - £c) 2 + 0(81) , 


(B- 5) 
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which in turn defines the inverse transformation, 

£ = £c± V2CRc(M 2 c - 1 )' 1/4 \/FSK • ( S -6) 

This implies that a point S in the neighborhood of a branch J\ (B.l), on the side with 

(7 c(^-Cc)>0, {B.l) 

is mapped to two points on both sides of J . For a point 5 on the other side of J , with 
(£ _ £ c )cr' < 0, (B.6) says only that there is no inverse solution S in the neighborhood of 
point S c - Because of (B.l), (B.6) does not include inverse solution(s) of (1.12) which are 
away from the branch J\. From (2.6), <r' c and are of oppsite sign and (B.6) implies 

(i) Along a branch of j where cr c increases (decreases) with increasing r c , a point 5 on 
either side of J is mapped to a point 5 on the right (left) hand side of J with £ - £ c 
positive (negative) while £ c decreases (increases) along J. The multiplicity of the 
point S is at least 2. 

Now consider the case that the point 5(£, a) is close to an extremum point on J , say 
at point A'i (£i , u\) where a[ — 0. Here we use subscript 1 to denote the value at. rj. The 
distance between S and 7?i is small, i. e., \SKi \ = \/(£ - 6) 2 + (^ - ^i) 2 < 1, and hence 

C\T-n\ = \R-Ri\ < \SKx\ <1 . {B. 8) 

Because of o\ — 0, we need to include at least the third order terms in the representation 
of the transformation (B.3) by a power series of r - t x or rather R- R\, 

(-( I =i-( 1 -M,(R-R l )+^M[(,R-R 1 ) 1 -^M'^R-R 1 f+0((R-R 1 ) t ). (B. 9) 


With £ = VR 2 — cr 2 , we expand £ — £1 i n a double power series of R — i?i and a — <j\. 
Using (2.4b), we note 

|<r - cri | < |<7 C - crj | = 0{{ic ~ 6) 2 ) > (5.10) 

while \R-Ri \ can be 0(|£ c - £ 1 1). We keep 0{\R - Ri | 3 ) terms but not fourth order terms, 
e. g., {a - ai) 2 and (R - R\f{o - <7i), and get 


f - 6 = (|m* - «.) + - x>) 2 + \(§mr - R>f 

+ -<n) + (J^h(R - RiK " - ffi) + 0((R -Ri) 1 ). 


With the partial derivatives at t\, 

gi = d = Mi w- 1 ) 

dR ~ & 1 ’ ai ? 2 
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_ _ZL 

da Ci 
(B.ll) and (B.9) yield 


5 3 C _ 3 a\R x _ 3M 3 (M 2 - 1) 

ai ? 3 g £ 2 



and 


g 2 C 

dadR 


a x Ri _ MjVM 2 - 1 
8 “ 


C - Cl = - \/^l - 1 (<T - <7i) + 
<M 2 y/M 2 - 1 


A/ 2 \JM\ - 1 


(<7 - <Ti )(R ~ #l) 


(£. 12 ) 

6C^— («-^) s + 0(Kc-Ci| 4 ). 

Since <7, <7i and i?i are known, (B.12) becomes a cubic equation for (i? — i?i) or R, 

(R - Ri) 3 + P(R - Ri) + g = 0 , (£.13) 


where 


P = 


— 6C 2 (a — a i) 


Q = 


6C 2 RiH 


a'l M 2 ^ M 2 - 1 


and H = (t-Z 1 ) + (v-<r 1 )y/Mjj -1 . 


Here £ denotes the horizontal distance from point (C, <r) to the tangent line of J at the 
cusp point K\. The character of the roots of (B.13) depends on the discriminant of the 
cubic, 0 = q 2 / 4 + p 3 / 27. If ft < 0, (B.13) has three distinct real roots, 0 = 0, three real 
roots with two equal and 0 > 0, only one real root. 

In order to have three real roots it is necessary for p < 0 that 


(a — ai)a'l > 0 . 


(*.14) 


This says that 

(ii) For the cubic (B.13) to have three real roots, it is necessary that a should be below 
the local maximum <Ji on J or above the local minimum. 

When ft > 0, we have 


2 „ SC*Mi(M? - 1 ) , 

H ~ (< ’-‘ n) 


(*•15) 


With the equality sign, (B.15) becomes the equation describing the two branches of J near 
the cusp S\. For a point 5(C, a) in the interior bounded by these two branches, we have 
ft > 0, and hence there are three distinct real roots for C or r. In the exterior there is only 
one real root. Thus using (B.14) and (B.15) we conclude that 

(iii) For a cusp with a local maximum a\, a point S in the interior bounded by these 
two branches has image point S inside the domain T> + bounded by J and two image 
points in V~ lying to the right and left of the local maximum point Si . For a point 
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S exterior to the region bounded by the two branches, there is only one image point 
Sint)-. 

(iv) For a cusp with a local minimum <72, a point S in the interior bounded by the two 
branches has one image point S in T>~ and two image points in Z>+ lying to the right 
and left of the local minimum point S 2 - For a point S exterior to the region bounded 
by the two branches, there is only one image point S in T> + . 

Again we note that the above analysis is local and does not account for the roots of (1.12) 
not in the neighborhood of an extremum point, if any. Thus we say 

(v) The multiplicity of a point 5(£, a) near a cusp and in the interior bounded by the two 
branches is at least 3. 
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Fig. 1 The coordinate system (x, y, z) with the medium at rest and the system ( x , y, z), with 
x = x + X(t). 



Fig. 2 The domain of dependence T> of point P at instant t > 0. in a meridian plane, the £cr 
plane, with origin P . With length scaled by Cf , the domain T) is the unit semicircular 
disk above the £axis. A fixed point S shown at t = 1 is moving towards the origin 
as t increases. Also shown is the partition of T) into two nonoverlapping regions, 
by the radial line PT at ^ = arccos 1/M representing the critical curve for a constant 
supersonic speed at Mach number M = \/5. 
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Fig. 3 The domain V and point P are the images of V and P in Fig. 2. The domain V, 
which depends on the relative velocity is shown here in the meridian £ a plane 
with origin P for a constant supersonic speed at M = with length scaled by Ct. 


f/(Ct) 

r 



Fig. 4 The partition of the domain of dependence V of point P in a meridian £<7 plane for 
a constantly accelerating motion. Here U(t) = Ct , r € [0, y/E], with final Mach 
number M = y/b. The critical curve J has a maximum ct(ti) at point K\ with 
ti = 5 1,/6 . The point K\ divides J into J\ and The critical curve partitions f> 
into V±. The circular arc C+ partitions V + into Fj -1 " and V^~ . 
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Fig. 5 Images of curves shown in Fig. 4 in the £<7 plane. They are the semicircle B, two 
branches of the critical curve J, J\ and J 2 , and the circular arc Cf through the cusp 
A'i on J . An enlargement of the neighborhood of the cusp K\ is shown in the insert. 
Also shown are the common tangent to the two branches, and the tangents at points 
on the branches. These tangent lines are the local Mach lines. 
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Fig. 6 The partition of the domain of dependence f> of point P in a meridian £<j plane for a 
constant rate acceleration from a low supersonic speed. Here M(r) = 1.15 + r 2 , r £ 
[0,2.25], with final Mach number M — 6.2125. The critical curve J has a local 
maximum at point K\ and minimum at point /\ 2 - The curve partitions V into V ± . 
y + is partitioned by the circular arc C+ starting from point I\\ into V r 1 + and F 3 + , and 
V~ by arc from point /? 2 into and V± . 
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Fig. 7 Images of curves shown in Fig. 6 in the £<7 plane. They are the semicircle 3. the 
three branches of the critical curve, Ji, J 2 and J3, and two circular nrcv 1’* and 
from the two cusps K\ and I \ 2 on J . Note the small difference between C 2 . and 
the semicircle B. A 20 x enlargement of the neighborhood of the cusps J\ 2 and A'3 
is shown in the insert. Note that the branch Jz begins at point T tangent to the 
semicircle 3 and continues downward to the right hand side of 3. 
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Fig. 9 The partition of the domain of dependence t> of point P in a meridian fcr plane for 
an accelerating, decelerating and accelerating supersonic speed starting from rest at 
t = 0. The subsonic and supersonic intervals are [0, 1) and (1,5] with the deceleration 
in (1.7,3.05). The critical curve J has a local maximum at point Ki , a minimum at 
I< 2 and then a maximum at K 3 . The curve J partitions t> into V ± . V~ is partitioned 
by the circular arc C 2 starting from point K 2 into two subdomains V 2 ~ and , and 
V + by arcs from point K\ and from point K 3 into three domains V)" 1 ", V' 3 + and 

. 


35 



Fig. 10 Images of the five nonoverlapping subdomains of V, shown in Fig. 9. are superimposed 
onto the £cr plane in (10a). Shown in Fig. (10b) is the 5x enlargement of the region 
containing the four branches of the critical curve, J it i - 1,2, 3, 4. The numerals in 
bold face indicate the multiplicities of the regions bounded by 3i-,B and/or the /-axis. 







Fig. 12 Images of the four nonoverlapping subdomains of f> shown in Fig. 11 are superposed 
on to the plane in Fig. 12a. Shown in Fig. 12b is the 45 x enlargement of the region 
containing the two cusps. 












Fig. 15 The partition of the domain of dependence T) of point P at time t with two supersonic 
intervals, M(r ) = Mo — br + t 2 / 2 for 0 < r < t, with Mo and t the same as in Fig. 13 
but b = 0.8. 
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Fig. 19 Propagation of signals from point S initiated at different 
vation point P at time t, as explained in Appendix A. 


retarded times to the obser- 
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